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ABSTRACT 

We investigate the emission properties of polycyclic aromatic hydrocarbons (PAHs) in various metal- 
hcity environments with the Infrared Spectrograph on board Spitzer. Local giant H ii regions are used 
as references as they enable access to the distinct interstellar medium components that contribute to 
the mid-infrared spectrum of star-forming galaxies: photodissociation regions (PDRs), photoionized 
gas, stellar clusters, and embedded regions. Three objects are considered, NGC 3603 in the Milky 
Way, 30Doradus in the Large Magellanic Cloud, and N66 in the Small Magellanic Cloud. From the 
variations of the PAH/14 /xm ratio, we find that PAHs are destroyed in the ionized gas for a radiation 
field such that [Ne iii]/[Ne ii]> 3. From the variations of the PAH/Hua ratio, we find that the PAH 
emission sources in the giant H II regions follow the same photodestruction law regardless of nietal- 
licity. We then compare these results with observations of starburst galaxies, H ii galaxies, and blue 
compact dwarf galaxies (BCDs). While the integrated mid-infrared spectra of BCDs are reminiscent 
of a warm dusty ionized gas, we observe a significant contribution to the PAH emission in starburst 
galaxies that is not arising from PDRs. 

Subject headings: HII regions, ISM: dust, ISM: atoms, ISM: molecules, infrared: ISM, Telescopes: 
Spitzer, (ISM:) supernova remnants, galaxies: star clusters 



I 1. INTRODUCTION 

From a simplistic chemical viewpoint, the abundance 
of molecules in the interstellar medium (ISM) of galax- 
ies should scale with the metallicity of the environ- 
ment, either because of the presence of nucleosynthetic 

[ constituents or because of the need for catalysts (dust 
grains). Obviously, indirect consequences of the low 
metal abundance also impact molecule survival, such as 

' the lack of certain coolants (such as [O i] or CO) that 
results in a globally warmer ISM, and hard interstellar 
radiation field (ISRF) due to energetic photons arising 
from low-metallicity ionizing stars. 

Observations show that nearby star-forming galaxies 
with metaUicities lower than ss 1/6 (12-hlog (O /H) < 
7.9) seem to lack CO molecules in giant molecular clouds 

', (GMCs; Arnauh et al. 1988; Taylor et al. 1998) and poly- 
cyclic aromatic hydrocarbons (PAHs) in photodissocia- 
tion regions (PDRs) (Engelbracht et al. 2005; Madden et 
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al. 2006; Wu et al. 2006; O'Halloran et al. 2006; Rosen- 
berg et al. 2006). Even H2, although it is detected in 
PDRs of metal-poor star-forming galaxies (Vanzi et al. 
2000; Hunt et al. 2010), is not seen in absorption in their 
diffuse ISM (Vidal-Madjar et al. 2000; Lebouteiller et al. 
2004; Thuan et al. 2002, 2005). This resuh iUustrates the 
indirect effect of low-metallicity since the paucity of dif- 
fuse H2 is related to the combination of low H i volumic 
density, dust grain abundance, and hard UV radiation 
(Vidal-Madjar et al. 2000). 

In the present paper, we study the emission proper- 
ties of PAHs, which are carbon-based molecules sug- 
gested by Leger & Puget (1984) and AUamandolla et 
al. (1985) to explain the mid-infrared (MIR) emission 
bands seen ubiquitously toward star-forming regions. 
PAH features are often used as a diagnostic tool, as it 
is thought that PAH emission bands could trace star- 
formation on timescales on the order of 0-stars lifetime 
(Forster Schreiber et al. 2004; Weedman & Houck 2008) 
or B-stars lifetime (Peeters et al. 2004). It is therefore 
essential to understand the mechanisms of their forma- 
tion and destruction as well as the physical conditions in 
which PAH emission occurs. 

The weakness of PAH emission in low-metallicity star- 
forming galaxies is thought to be due to a combination 
of several effects: insufficient chemical evolution of the 
environment leading to the lack of PAH progenitors (hy- 
drocarbons) , enhanced molecule destruction by energetic 
photons and by shocks, and smaller volume of PAH- 
dominatcd regions. The main difficulty in investigating 
these parameters is their cross-correlation. Although a 
trend of decreasing PAH emission with both metallicity 
and ionizing radiation hardness has been found in sev- 
eral studies (Madden 2000; Madden et al. 2006; Wu et 
al. 2006; Engelbracht et al. 2006), the dominant param- 
eter remains unknown. The study of H ii regions within 
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MlOl by Gordon et al. (2008) suggests that the aro- 
matic band emission tends to depend more on the ioniza- 
tion of the environment rather than on metallicity. The 
photodestruction mechanism could be isolated by inves- 
tigating PAH emission at ionization fronts in reflection 
nebulae (e.g.; Sellgren et al. 1985) and across spatially- 
resolved H II regions where the metallicity is known to 
be uniform (Madden et al. 2006; Lebouteiller et al. 2007; 
Bernard-Salas et al. 2011). 

The present study aims at better understanding how 
the PAH emission in galaxies depends on the ISM 
morphology, by comparing the MIR properties of star- 
forming galaxies with those of giant H ii regions. The 
latter enable access to the distinct physical regions con- 
tributing to the composite mid-infrared (MIR) spectrum 
of star-forming galaxies, such as warm photoionized gas, 
PDRs, molecular clouds, shock fronts, stellar clusters, or 
supernovse remnants. Three well studied giant H ii re- 
gions with different metallicities are used as templates. 
The objects are NGG 3603 in the Milky Way with essen- 
tially a solar metallicity, 30Doradus (hereafter 30 Dor) 
in the Large Magellanic Cloud (LMC) with 0.6 Zq, and 
N66/NGG346 in the SmaU Magellanic Cloud (SMC) 
with 0.2 Zq. Although these regions might not be repre- 
sentative of typical H ii regions in star-forming galaxies, 
it is possible to distinguish their ISM components in the 
MIR to compare to integrated spectra of galaxies. 

This paper is part of a series meant to investigate the 
MIR properties of nearby giant H ii regions. Lebouteiller 
et al. (2007) analyzed small spatial-scale variations of the 
MIR features in NGC 3603 to characterize PAH destruc- 
tion by energetic photons. A valuable advantage is that 
chemical abundances are remarkably uniform across the 
region, hypothesis confirmed a posteriori by Lebouteiller 
et al. (2008). Hence, it provided a unique possibility to 
isolate the metallicity parameter and study the specific 
effect of ionizing radiation hardness. The two other ob- 
jects are being analyzed a similar way (Bernard-Salas et 
al. 2011; Whelan et al. 2011). 

We focus here on the brightest regions seen in the 
MIR to compare the PAH properties in the 3 giant H ii 
regions and to understand how PAH emission depends 
on both the physical and chemical conditions (Sect. U]). 
We then view our results alongside integrated spectra of 
star-forming galaxies to determine general diagnostics of 
molecular and dust properties in various environments 
(Sect.[S|). The two specific issues we wish to address con- 
cern the infiuence of metallicity on the PAH-to-dust ra- 
tio, and whether the relative contribution of ionized gas 
versus PDRs in the integrated spectra of star-forming 
galaxies can explain the observed PAH deficiency in low- 
metallicity objects. Finally, we examine in Section [5] the 
relation between PAH emission and molecular hydrogen. 

2. PAH EMISSION PROPERTIES 

PAHs represent an important component in dust mod- 
els (e.g.; Draine & Li 2001) and their importance is mul- 
tifold (e.g., Tielens et al. 2008; GaUiano et al. 2008b). 
Most notably, they allow efficient cooling of the ISM, 
while providing most of the heating in PDRs (Bakes & 
Tielens 1994; HoUenbach & Tielens 1997). In addition, 
they contain a significant fraction (< 20%) of the carbon 
depleted from the gaseous phase (Allamandola, Tielens 
& Barker 1989; Dwek et al. 1997; Zubko et al. 2004). We 



briefiy review here the parameters responsible for the 
production and destruction of PAH molecules. We also 
emphasize the importance of the ISM morphology on the 
detection of PAH emission. 

2.1. Formation 

Carbon-rich asymptotic giant branch (AGB) stars pro- 
vide a priori all the necessary ingredients (high collisional 
rates, cooling flow, shielding, periodic shocks) to produce 
significant amounts of PAH progenitors (Latter 1991). It 
is now thought that aliphatic hydrocarbons, which have 
been detected in large amounts in the envelopes of AGB 
stars (Cernicharo et al. 2001; Sloan et al. 2005; Boersma 
et al. 2006; Bernard-Salas et al. 2006), can lead to the for- 
mation of aromatic molecules via photoprocessing of hy- 
drocarbon bonds. Theoretical studies support this PAH 
formation mode (e.g., Cadwell et al. 1994; Cherchneff 
2006; Joblin et al. 2008). 

With this scenario in mind, the PAH formation rate 
over a galaxy history could be mostly a function of the 
rate of low-mass stars entering the post- AGB phase. The 
spectral energy distribution modeling of galaxies by Gal- 
liano et al. (2008b) shows that PAH production in AGB 
stars is consistent with the PAH abundance evolution 
with metallicity (see also Mufioz-Mateos et al. 2009). 
In particular, a young galaxy experiencing its first star- 
formation episode should contain little or no PAHs up 
until w 500 Myr. Based on this argument, Galliano et 
al. (2008b) proposed that delayed injection of PAHs by 
AGB stars is responsible for the low PAH abundance in 
metal-poor galaxies. If the dominant production process 
is indeed due to AGB stars, the PAH formation rate de- 
pends indirectly on the observed ISM metallicity, and the 
PAH abundance is expected to peak around solar metal- 
licity sources (Galliano et al. 2008b), marking the tran- 
sition from carbon-rich to oxygen-rich AGB stars (e.g., 
Renzini & Voh 1981). 

It must be noted that several other mechanisms have 
been proposed for PAH formation in the ISM, most no- 
tably splitting of graphitic grains in shocks (Tielens et 
al. 1987; Jones, Tielens & HoUenbach 1996), accretion 
of C+ in the diffuse ISM (Omont 1986, Puget & Leger 
1989), gas- formation by ion- molecule reaction in dense 
clouds (Herbst 1991), or photoprocessing of interstellar 
dust mantles (Greenberg et al. 2000). Most of these pro- 
cesses unfortunately still lack observational constraints. 

2.2. Destruction mechanisms 

Molecule destruction eventually controls the PAH 
abundance in a galaxy, especially in extreme environ- 
ments (GaUiano et al. 2003; 2005; Madden et al. 2006). 
Metal-poor starburst galaxies globally exhibit harder ra- 
diation fields than more metal-rich galaxies since the ra- 
diation from ionizing stars have relatively lower opac- 
ity (e.g., Schaller et al. 1992; Schaerer & Maeder 1992; 
Schaerer et al. 1993; Charbonnel et al. 1993). As a result, 
the ISM is relatively warmer and the possible deficiency 
of some coolants (in particular CO) could exacerbate 
this effect. Furthermore, if massive stars in metal-rich 
starburst galaxies evolve in heavily extinguished regions 
(even opaque to MIR wavelengths), the observed inter- 
stellar radiation field (ISRF) would be found to be softer 
than in metal-poor galaxies (Rigby & Rieke 2004). 
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Radiation in low-metaflicity environments is poten- 
tially energetic enough and intense enough for photons 
to destroy PAH molecules, depending on the molecule 
size (Omont 1986). PAHs can be heated to tempera- 
tures of > 1 000 K by a single photon event. The energy 
deposited and thermalized into the molecule can lead to 
evaporation of bonds. Direct photolysis by UV photons 
is also an important process in the vicinity of young mas- 
sive stars. 

Observationally, PAH emission is seen closer to the ion- 
izing source compared to H2 and CO (e.g., Hollenbach 
& Tielcns 1999; Cesarsky et al. 1996). Photodestruc- 
tion of PAHs effectively occurs in the ionized gas, where 
the observed intensity of PAH bands depends strongly 
on the physical conditions of the gas (e.g., Giard et al. 
1994; Kassis et al. 2006). As opposed to H2 and CO, 
PAH molecules are not self-shielded which, depending 
on the physical conditions, can lead to full inhibition of 
PAH emission at PDRs while H2 and CO are still able 
to emit. 

Finally, it must be added that sputtering can occur 
in strong shocks (e.g., Micelotta, Jones & Tielens 2009) 
and in the hot gas, providing yet another means to de- 
stroy PAH molecules. However, shocks seem to be able 
to destroy all the dust content (Reach et al. 2002), so 
that they might not be good candidates to explain the 
observed deficiency of the PAH emission relative to dust 
continuum emission. In addition, rather than destroying 
PAHs, shocks could in fact result in molecule formation 
due to the splitting of large graphitic grains fSect. l2.11 ). 

2.3. PAH emission 

PAH emission arises in the neutral gas of galaxies, ei- 
ther diffuse (Haas et al. 2002) or dense (i.e., in PDRs, 
see introduction). Observations of integrated galaxies are 
more likely to show a significant contribution from PAHs 
in the diffuse ISM as compared to single star-forming 
region, where PAH emission mostly comes from PDRs, 
thus depending on the distribution of molecular clouds 
in the vicinity of ionizing stars. 

Both the PAH emission and the warm dust emission 
trace on first order the presence of young massive stars. 
However, while warm dust emits in the ionized gas of 
H II regions (the grains in thermal equilibrium becomes 
hotter as the energy density U increases. Sect. 14. 4p . PAH 
emission is intense at the vicinity of dense clouds such as 
GMCs. This different spatial origin can have important 
consequences in dwarf galaxies since these are prone to 
a lack of GMCs. In the SMC, Leroy et al. (2007) found 
that GMCs are either small or translucent with an av- 
erage extinction several magnitudes lower than clouds in 
the Milky Way. Because of the low volume filling fac- 
tor of GMCs combined with the beam size, the CMC 
properties in more distant objects is much less known. 
GMCs have been detected in nearby starburst galaxies 
and a few dwarf irregulars while non-detections are fre- 
quent in blue compact dwarf galaxies (BCDs) because of 
the low CO abundance in metal-poor environments that 
severely reduces the self-shielding of the CO core (e.g., 
Gondhalekar et al. 1998 and references therein). 

Although dwarfs contain H i clumps with similar sizes 
and masses as the GMCs in the Milky Way (Lo et al. 
1993), most of their H i mass is expected to lie in the 
warm neutral medium. Indeed, in galaxies with a high 



star-formation rate, the low-density hot bubbles created 
by stellar clusters through stellar winds and successive 
SNe could well overlap so that the cold neutral medium 
is finally disrupted by the stellar feedback. PAH emis- 
sion should thus be hampered since PAH molecules are 
destroyed in the warm ionized gas and in the hot gas 
(Sect.[13. 

As a conclusion, the PAH abundance, when compared 
to the dust component (warm or cold), strongly depends 
on the ISM morphology and in particular the distribu- 
tion and number of molecular clouds. In peculiar ob- 
jects such as BCDs, the observed lack of PAHs might 
thus not be a strong constraint on their actual forma- 
tion/destruction over the galaxy history. Because of the 
dependence on the ISM morphology, only comprehensive 
galaxy chemical evolution models including GMC forma- 
tion and stellar feedback would allow us to understand 
how the galaxy mass and the starburst properties might 
affect the integrated PAH emission in galaxies. In this 
study, we address this issue empirically by comparing the 
emission properties of single massive H 11 regions with in- 
tegrated star-forming galaxies. 

3. DATA ANALYSIS 
3.1. Observations 

The observations of the giant H 11 regions are presented 
in Lebouteiller et al. (2007; 2008) and we summarize here 
the most important aspects. NGC3603, 30 Dor, and 
N66 were observed by the Infrared Spectrograph (IRS; 
Houck et al. 2004) onboard the Spitzer Space Telescope 
(Werner et al. 2004) with the high- and low-resolution 
modules. This resulted in spectral coverage from 5 fim 
to 35 fim. Data were processed using the Spitzer Science 
Center SI 3. 2 pipeline. We used the "basic calibrated 
data" product. Rogue pixels and on-the-fly flagged data 
were removed using IRSCLEAINH. Spectra were then ex- 
tracted using the SMARTO package (Higdon et al. 2004; 
Lebouteiller et al. 2010). 

In contrast with Lebouteiller et al. (2007) and Bernard- 
Salas et al. (2011) who used the spatial information 
available along the low-resolution aperture, we consider 
in this study the spectra of the nominal positions for 
which both high- and low-resolution modules are avail- 
able. The extraction aperture is 4.7" x 11.3" (full slit) for 
the Short-High module, 11.1" x 22.3" (full slit) for Long- 
High, 7.2" X 3.7" for Short-Low, and 20.4" x 10.7" for 
Long-Low. Spectra were normalized to the Short-High 
spectrum based on matching the dust continuum and as- 
suming a uniform light distribution within the apertures. 
Final spectra can be found in Lebouteiller et al. (2008). 
Since we consider only the nominal positions, there is 
an observational bias toward the MIR brightest sources 
within each giant H 11 region. There are no dedicated 
observations of quiescent regions, except the positions 
N66#l, #2, and #3 which are mispointings. 

The IRS spectra of giant H 11 regions probe vari- 
ous stellar/ISM components (Lebouteiller et al. 2007; 
Bernard-Salas et al. 2011), with most notably: 
Central star clusters. A strong continuum has been 

The IRSCLEAN software can be downloaded from 
http://ssc. spitzer. caltech. edu 

11 The SMART package can be found at 
http. ■/ /is c. astro. Cornell. edu/Smart/ 
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detected shortward of 8 fim in NGC 3603. The MIR spec- 
trum is otherwise similar to that of the warm photoion- 
ized gas. Spitzer/IRAC images show a dust- free cav- 
ity around the central clusters of NGC 3603 and 30 Dor, 
probably carved out by the winds of massive stars, with 
little or no MIR emission. Silicate dust emission was de- 
tected in emission toward, or close to, the central clusters 
of the 3 giant H ii regions, hinting at low optical depth 
lines of sight (Whelan et al. 2010). 

Warm photoionized gas. It is detected via the H i 
recombination line Hua and the forbidden lines [S ill], 
[S iv] , [Ne ii] , [Ne iii] , [Ar iii] , and [Fe ill] . The ionized 
gas component harbors a strong and steeply increasing 
dust continuum longward of 10 /xm. 
PDRs. Each giant H li region contains several indi- 
vidual PDRs. These are located further away than the 
ionized gas counterpart from the central cluster, indicat- 
ing that the radiation coming from the central cluster 
dominates the energetics. PDRs show bright PAH bands 
as well as significant H2, [Ar 11], [Si 11], and [Fc 11] line 
emission. 

Embedded MIR bright regions. Embedded source 
show silicate dust in absorption. Only 30 Dor and N66 
show sources with deep silicate absorption in our obser- 
vation. PAHs are also detected toward these regions in 
some cases. 

Since we are dealing with extended sources, the extrac- 
tion aperture (Sect. [33} probes larger physical regions 
in more distant objects. Therefore, while observations 
of NGC 3603 (« 7kpc away, e.g., Moffat et al. 1983) 
can be centered on a single PDR, observations of 30 Dor 
(« 49kpc, e.g., Selman et al. 1999) and N66 (w 61kpc, 
e.g., Hilditch et al. 2005) can include several PDRs. As 
an illustration, the projected size of the IRS Short-High 
aperture is 0.4 pc in NGC 3603, 2.7 pc in 30 Dor, and 
3.3 pc in N66. Furthermore, projection effects prevent 
us from asserting with accuracy the physical association 
between the tracers of the ionized gas and those of the 
PDRs. Nevertheless, observations were designed to study 
the brightest MIR knots which are dominated by a single 
region within the aperture with little contribution from 
nearby regions (see also Lebouteiller et al. 2007). 

In the following we compare the giant H 11 regions with 
a sample of BCDs (Wu et al. 2006), starburst galaxies 
(Brandl et al. 2006; Bernard-Salas et al. 2009), and H 11 
galaxies (Smith et al. 2007). 

3.2. Tracers 
3.2.1. Molecules and dust 

We measure PAH emission by summing the flux of the 
most prominent bands in the IRS wavelength range, i.e., 
5.7 /im, 6.2 /im, 7.7 /im, 8.6 /im, and 11.3 /im. PAH fea- 
tures were measured using the PAHFIT profile decompo- 
sition algorithm (Smith et al. 2007). The 11.3 /xm PAH 
flux was taken as the integrated fiux of the 11.23 /im 
and 11.33 /im PAHFIT features. Similarly, the 7.7 /im 
PAH fiux represents the integrated flux of the 7.60 /im 
and 7.85 /im PAHFIT features. The H i recombination 
lines at 7.46 /im (Pfa) and 7.50 /im (Hu/3) contaminate 
the 7.7 fim PAH complex by up to 5%. We added these 
lines manually in PAHFIT to obtain a satisfactory fit of 
the PAH complex. The plateau underneath the 7.60 /im 
and 7.85 /im features was included in the fit but ignored 



in the calculations of the PAH band fiuxes, as it prob- 
ably originates from larger molecules such as PAH clus- 
ters with several 100s carbon atoms (e.g., Bregman et al. 
1989). PAH fiuxes for the giant H 11 regions are reported 
in Tabled! 

While Smith et al. (2007) measured PAH emission with 
PAHFIT in H 11 galaxies, a spline method was originally 
used to remove the underlying spectral continuum in the 
BCDs and in the starbursts. Hence for consistency, we 
re-calculated PAH fluxes in these objects using PAHFIT 
(Table[2]). As expected, the new PAH fluxes are different 
than those measured by Wu et al. (2006), Brandl et al. 
(2006), and Bernard-Salas et al. (2009), while variations 
in band ratios are mostly unchanged. A discussion on the 
comparison between the spline method and the profile 
decomposition method can be found in Galliano et al. 
(2008a). 

In the following, we compare the PAH emission to that 
of the dust grains responsible for the continuum emis- 
sion at wavelengths longward of 10 /im. The grain emis- 
sion is determined by integrating the continuum flux in 
a 1 /im wide window centered at 14.3 /im, where there 
is no contamination by lines or molecular features. The 
continuum flux is given in Table[TJ Close to the ion- 
izing sources, the continuum emission is expected to 
be dominated by warm big grains (> lOOA radius) m 
thermal equilibrium with the ISRF. At larger distances, 
i.e., at PDRs, dust particles are colder and there might 
be a contribution of smaller grains (lO-lOOA) stochas- 
tically heated (Desert et al. 1990). Note however that 
this paradigm is modulated by the ISRF intensity and 
electron density that control the heating mechanism (see 
Sect.Ol). In NGC 3603, Lebouteiller et al. (2007) ob- 
served that, unlike PAHs, warm grains traced by the 
14 fim continuum are able to survive the radiation fleld in 
the ionized gas. The grain emission in this region is spa- 
tially correlated with [S iv] emission, little with [Ne 11], 
and is anticorrelated with [Ar li]. Accordingly, warm 
grains are the main dust component in the ionized gas 
(see also Cesarsky et al. 1996; Verstraete et al. 1996). 
The same result is found by Bernard-Salas et al. (2011) 
in 30 Dor. 

Finally, we use the H2 S(l) and S(2) lines, at respec- 
tively 17.03 /im and 12.28 /im, to trace molecular hydro- 
gen (Tablelll). The S(0) line at 28.22 /im is not detected 
toward most positions. In H 11 nuclei, the S(0) line con- 
tributes to about 25% of the total of S(0-2) (Roussel et 
al. 2007). 

3.2.2. Ionized gas 

Several forbidden emission-lines can be used to trace 
the warm photoionized gas and its properties (ISRF in- 
tensity and hardness). On one hand, the H i recombina- 
tion line 12.37 /im Hua depends on the physical condi- 
tions (electronic density and temperature) but not on the 
ISRF hardness or on the metallicity. On the other hand, 
the sum of the 12.81 /im [Ne 11] and 15.56 /im [Ne ill] line 
fluxes depends on the neon abundance but not on the 
radiation hardness as Ne 11 and Ne ill are the dominant 
ionization stages of neon in the warm ionized gas. The 
quantity we use in the following as a proxy for the ex- 
cited neon flux is [Ne 11] -I- C [Ne ill] where C accounts 
for the fact that the 2 lines have different spontaneous 
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TABLE 1 

PAH AND CONTINUUM FLUXES IN GIANT H II REGIONS. 
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Position 


PAH 5.7 


"r> ATT r' o . 

PAH D.2 /jm 


Ti A TJ ^7 ^ . 

PAH 7.7 /im 


T~> A TT O 1^ . 

PAH 8.6 


T^ A T T 1 "1 O . 

PAH 11.3 


/(14^m) 


NGC3603#1 


0.22 (0.51) 


12.57 (0.91) 


20.18 (0.51) 


8.6 (0.51) 


11.16 (0.78) 


22.01 


NGC3603#3 


43.01 (0.83) 


301.09 (0.95) 


678.82 (0.83) 


210.26 (0.83) 


401.7 (3.5) 


799.55 


NGC3603#4 


51.53 (20.12) 


320.55 (15.62) 


500.42 (20.12) 


159.3 (20.12) 


345.09 (19.53) 


1332.09 


NGC3603#5 


52.71 (0.85) 


418.98 (0.95) 


947.76 (0.85) 


289.61 (0.85) 


323.87 (3.33) 


865.57 


NGC3603#6 




15.32 (1.2) 




11.88 (0.76) 


14.39 (3.79) 


309.55 


NGC3603#7 


12.15 (1.94) 


35.89 (1.4) 


24.21 (1.94) 


35.24 (1.94) 


8 (0.68) 


8371.07 


30Dor#2 




10.03 (0.31) 


25.17 (0.34) 


8.86 (0.79) 


15.3 (0.47) 


10.27 


30Dor#3 


1.13 (0.63) 


4.15 (0.44) 


5.08 (0.63) 


1.94 (0.63) 


4.61 (0.6) 


88.02 


30Dor#4 


2.52 (0.47) 


3.62 (0.36) 


2.37 (0.47) 




9.85 (0.22) 


230.33 


30Dor#5 


4.34 (0.51) 


38.11 (0.41) 


84.69 (0.51) 


29.71 (0.51) 


40.45 (0.76) 


33.74 


30Dor#6 


3.27 (0.77) 


40.56 (0.65) 


83.25 (0.77) 


26.62 (0.77) 


58.58 (0.92) 


156.98 


30Dor#7 


1.43 (0.59) 


16.87 (0.54) 


23.78 (0.59) 


9.23 (0.59) 


20.05 (0.65) 


308.08 


30Dor#8 


9.56 (0.64) 


67.53 (0.54) 


153.54 (0.64) 


46.96 (0.64) 


85.06 (0.86) 


224.46 


30Dor#10 


8.38 (0.58) 


63.08 (0.47) 


145.14 (0.58) 


45.27 (0.58) 


70.42 (0.75) 


133.5 


30Dor#ll 


3.52 (0.68) 


26.44 (0.4) 


48 (0.68) 


19.84 (0.68) 


41.64 (0.79) 


168.71 


30Dor#12 


3.23 (0.62) 


30.26 (0.53) 


75.88 (0.62) 


22.53 (0.62) 


43.29 (0.88) 


140.84 


30Dor#13 


2.25 (0.56) 


13.17 (0.38) 


9.7 (0.56) 


0.67 (0.56) 


6.15 (0.37) 


38.14 


30Dor#14 


4.5 (0.6) 


48.79 (0.51) 


120.79 (0.6) 


36.38 (0.6) 


79.95 (1) 


467.99 


30Dor#15 


3.59 (0.63) 


58.87 (0.46) 


124.11 (0.63) 


32.65 (0.63) 


87.99 (0.99) 


635.24 


30Dor#16 




16.39 (0.37) 


29.9 (0.43) 


7.44 (0.19) 


18.35 (0.82) 


466.53 


30Dor#17 


10.08 (1.11) 


30.04 (0.62) 


71.59 (1.11) 


52.68 (1.11) 


53.56 (1.02) 


176.05 


N66#l 




2.27 (0.12) 


4.55 (0.46) 


0.33 (0.30) 


3.39 (0.19) 


4.25 


N66#2 


0.18 (0.12) 


3.54 (0.09) 


1.33 (0.12) 


1.59 (0.12) 


5.48 (0.16) 


2.35 


N66#3 


0.53 (0.12) 


1.4 (0.09) 




0.43 (0.12) 


1.18 (0.1) 


1.61 


N66#5 




2.88 (0.48) 


9.24 (0.89) 


1.12 (0.42) 


4.72 (0.62) 


39.02 


N66#6 


0.9 (0.55) 


3.62 (0.35) 




4.05 (0.55) 


4.27 (0.3) 


31.84 


N66#7 




2.81 (0.29) 


4.15 (0.56) 


0.75 (0.52) 


3.1 (0.41) 


5.72 


N66#8 


0.03 (0.4) 


2.84 (0.36) 


5.26 (0.4) 


1.18 (0.4) 


4.12 (0.43) 


4.84 


N66#9 




5.28 (0.34) 


9.94 (0.97) 


2.63 (0.35) 


7.93 (0.5) 


12.76 


N66#10 


0.08 (0.47) 


2.23 (0.32) 


5.22 (0.47) 


0.68 (0.47) 


2.77 (0.51) 


7.04 


N66#ll 


0.46 (0.52) 


1.41 (0.31) 


0.64 (0.52) 


0.52 (0.52) 


2.28 (0.27) 


7.19 


N66#12 


0.5 (0.42) 


1.3 (0.31) 




0.23 (0.42) 


0.49 (0.23) 


5.43 


N66#13 




2.36 (0.34) 


1.07 (0.25) 


0.52 (0.31) 


2.3 (0.52) 


3.08 



Note. — Fluxes arc given in x lO^^'^ W cm^^. The 14 /im continuum flux is integrated on a 1 /^m wide window centered at 14.3 /^m. The 
uncertainty is given in parentheses. We assume an uncertainty of 10% for the 14 ^m continuum due to tire instrument flux calibration. PAH fluxes 
were measured using PAHFIT (Sect . [3. 2. H . 

^ Sources labels and coordinates are defined in Lebouteiller et al. (2008). 



transition rates and diff'erent level populations: 

„_ \NeIII] ^[Nell] 9[NeII] , . 

^ — -\ ■ -1 ■ : 

■^[Nell] ^[Afe///] 9[NeIII] 

where, for a given transition, A is the wavelength, A is the 
Einstein coefflcient, and g is the ratio of the population 
of the level from which the line originates to the total 
population of the ion. Results from Lebouteiller et al. 

(2008) for the giant H n regions and Bernard-Salas et al. 

(2009) for the starburst galaxies indicate that C k, 0.48 ± 
0.01, the deviation being dominated by slightly different 
electron densities and temperature. The quantity [Ne ii] 
+ C [Ne III] is directly proportional to the number of 
Ne H and Ne hi ions. 

Line intensities in the giant H ii regions are from 
Lebouteiller et al. (2008). The [Ne in] line could not 
be measured in positions #7 and #8 in NGC3603. For 
#7 we estimated the [Ne ill] flux using the correlation 
between [Ne iii]/[Ne ii] and [S iv]/[S ill] (Beirao et al. 
2006; Bernard-Salas et al. 2009). For position #8 we 
used an empirical relation between [S iv]/[Ne ii] and 
[Ne iii]/[Ne ii] ratios (this study and Groves et al. 2008). 
For the star-forming galaxies, we use the line measure- 
ments from Wu et al. (2006) for the BCDs, Bernard-Salas 
et al. (2009) for the starburst galaxies, and Smith et al. 
(2007) for H ii galaxies. H i line measurements for BCDs 



were provided by Yanling Wu (private communication). 

3.2.3. Radiation field 

The gas excitation is usually probed by line ratios 
of adjacent ionization states, such as [Ne iii]/[Ne ii] or 
[S iv]/[S III]. These ratios are proportional to the radia- 
tion hardness with also a — smaller — dependence on the 
ionization parameter U (e.g., Morisset et al. 2004). We 
ran a grid through the photoionization code CLOUDY 
(version cOS.OO; Ferland et al. 1998) to estimate the 
[Ne ili]/[Ne II] ratio in various conditions. A single star 
was chosen with a temperature T between 25 000 K and 
50 000 K, the ionization parameter U spanning from —4 
to —1, and the inner radius to the illuminated slab span- 
ning from 0.1 to Iparsec. Results are plotted in Fig.[TJ 
We flnd that the [Ne lll]/[Ne ll] ratio is proportional to 
T X g^jjfj jg ^]^yg jxiostly dependent on the star tem- 

perature hence the ISRF hardness. 

We favor hereafter the line ratio [Ne ill] 
15.56 /im/[Ne ii] 12.81 /im which is little affected 
by silicate dust extinction. Moreover, the two lines arise 
in the same module of the IRS so that flux scaling is not 
needed. Since the factor C is essentially the same for all 
objects fSect. lHTT^ . we ignore it for the neon line ratio. 

3.2.4. Chemical abundances 
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TABLE 2 

PAHFIT AND CONTINUUM FLUXES IN STAR-FORMING GALAXIES. 





•^mir'' 


7 b 

^opt 


PAH 5.7^111 


PAH 6.2^tm 


PAH 7.7/im 


PAH 8.6 Aim 


PAH 11.3 Atm 


/(14Atm) 


BCDs 


CG0752 






0.8 ( 0.21) 


4.36 (0.17) 


14.06 (0.53) 


2.1 (0.13) 


3.53 (0.26) 


6.39 


Haroll 


0.56 


0.18 


5.11 


( 0.25) 


9.65 (0.12) 


16.3 (0.49) 


2.65 (0.38) 


4.11 (0.17) 


90.24 


IIZw40 


0.31 


0.28 




( ■•• ) 


1.94 (0.12) 


3.61 (0.84) 


1.05 (0.44) 


2.05 (1.1) 


79.17 


Mark 1450 


0.69 


0.21 


0.14 


( 0.04) 


0.3 (0.02) 


0.34 (0.25) 


0.19 (0.03) 


0.19 (0.03) 


1.99 


Mark 1499 


0.55 


0.29 


0.13 


( 0.03) 


0.78 (0.03) 


1.91 (0.07) 


0.58 (0.02) 


0.82 (0.04) 


0.84 


NGC1140 


0.83 


0.46 


1.16 


( 0.22) 


6.41 (0.17) 


15.42 (0.65) 


4.28 (0.12) 


7.18 (0.36) 


8.65 


NGC 1569 


0.48 


0.3 


1.03 


( 0.26) 


8.77 (0.2) 


17.77 (1.11) 


4.94 (0.31) 


10.9 (0.39) 


78.45 


UGC 4274 






0.67 


( 0.14) 


4.81 (0.09) 


11.85 (0.77) 


3.4 (0.08) 


5.11 (0.27) 


4.66 


UM461 


0.19 


0.13 




( ■■• ) 










1.74 


Starbursts 


IC 342 


1.86 


1.55 


0(1 


3) 


264 (0.58) 


569 (4.06) 


189 (1.55) 


303.4 (0.88) 


343.78 


Mark 52 






0.93 


( 0-3) 


14.3 (0.3) 


35.5 (0.87) 


9.07 (0.21) 


12.46 (0.36) 


34.23 


NGC 520 






7.34 


( 1) 


155 (1.54) 


528 (8.19) 


163 (1.77) 


263 (6.99) 


93.31 


NGC 660 






47.8 


( 1-78) 


428 (2.07) 


1136 (15.44) 


371 (3.98) 


493 (4.86) 


214.03 


NGC 1097 






4.8 ( 


0.58) 


98.5 (0.88) 


282 (5) 


61.5 (3.66) 


129.6 (6.36) 


93.61 


NGC 1222 


1.07 




1.89 


( 0.22) 


30.8 (0.2) 


101.7 (1.33) 


22.9 (0.35) 


38.6 (0.76) 


51.01 


NGC 1365 


0.98 






( ■■■ ) 


42.3 (0.56) 


135.4 (2.54) 


36.6 (1.13) 


70.3 (1.5) 


208.44 


NGC 1614 


2.88 




10.8 


( 0.78) 


63.4 (1.07) 


242.4 (5.54) 


44.7 (3.42) 


63.7 (6.14) 


145.85 


NGC 2146 


1.82 




55.4 


( 0.66) 


578 (1.86) 


1677 (11.46) 


397 (3.74) 


549 (8.6) 


334.08 


NGC 2623 






2.41 


( 0-2) 


23.1 (0.31) 


104.5 (1.35) 


15.9 (0.4) 


24.3 (0.39) 


45.94 


NGC 3310 


1.32 


0.33 


4.33 


( 0.43) 


50.9 (0.36) 


117.8 (1.27) 


29.7 (0.25) 


40.3 (0.38) 


31.9 


NGC 3556 


0.6 




0.97 


( 0.33) 


28.3 (0.24) 


81.1 (1.49) 


22.1 (0.43) 


27.3 (0.43) 


22.28 


NGC 3628 


0.51 


0.81 


52.3 


( 1-26) 


375 (2.3) 


994 (9.02) 


279 (4.24) 


333 (4.39) 


89.14 


NGC 4088 


0.74 




2.65 


( 0.26) 


22.9 (0.39) 


60.5 (1.15) 


19 (0.34) 


25.5 (0.88) 


19.13 


NGC 4676 


1.17 




1.68 


( 0-11) 


10.3 (0.13) 


27.53 (0.8) 




9.06 (0.14) 


11.12 


NGC 4818 


2.09 




8.56 


( 0.23) 


73.6 (0.41) 


194.7 (1.7) 


47.7 (0.61) 


74 (1.23) 


120.63 


NGC 7714 


1.17 


0.4 


3.79 


( 0.19) 


35.3 (0.29) 


95.7 (1.38) 


27.2 (0.45) 


28.3 (0.44) 


68.57 



Note. — Fluxes arc given in x 10 W cm . The 14 //m continuum flux is integrated on a 1 //m window centered at 14.3 /j.m. The uncertainty is 
given in parentheses. We assume an uncertainty of 10% for the 14 //m continuum due to the instrument flux calibration. PAH fluxes were measured 
using PAHFIT ("Scct. [3?2Jl . 

^ MetalUcity estimate (in solar units) using the neon abundance derived using MIR fine-structure Unes by Wu et al. (2008) for the BCDs and 
Bernard-Salas et al. (2009) for the starburst galaxies. The solar reference we used is 12 + log(Ne/H)— 7.84 (Asplund et al. 2006). 
^ Mctallicity estimate (in solar units) from optical studies, Bergvall & Ostlin (2002) for Haro 11, Guseva et al. (2000) for IIZw40, Izotov et al. (1994) 
for Mark 1450, Shi et al. (2005) for Mark 1499, Izotov U Thuan (2004) for NGC 1140, Kobulnicky et al. (1997) for NGC 1569, Pilyugin et al. (2004) 
for IC342, Pastoriza et al. (1993) for NGC 3310, Engelbracht et al. (2005) for NGC 3628. The solar reference we used is 12 + log(Ne/H)^ 7.84 and 
12 + log(0/H)^ 8.66 (Asplund et al. 2006). 



^ 1 .00 







U=-l 
U=-2 






U = -2.5 - 






U = -3 : 






^^^^ U = -3.5 






U = -4 : 











3.5 4.0 4.5 5.0 

Teff (xlOHK) 

Fig. 1. — [Ne lll]/[Ne ll] predicted by CLOUDY as a function 
of the radiation hardness and of the ionization parameter (tracks). 
See text for model details. The gray zone indicates the range of 
line ratios observed in the objects of this study. 

The metallicity is derived from the neon abundance 
calculated using MIR fine-structure lines (Wu et al. 2008; 
Bernard-Salas et al. 2009). This choice was motivated by 
the fact that about half of the galaxies we consider do 
not have referenced abundance measurements from op- 
tical spectra. Moreover, IR abundance determinations 
are less affected by dust extinction. Our results remain 



unchanged by using the optically derived abundances 
fSect. [5T5|) . Metallicity values are listed in Tabled 

4. GIANT H II REGIONS 

4.1. Physical conditions 

The observations of the three giant H ii regions span a 
wide range of metallicity, while at the same time probing 
different physical conditions within each region (Fig. [2]). 
N 66 is characterized by a globally harder ISRF than the 
two other regions. Furthermore, although NGC 3603 and 
30 Dor are significantly more metal-rich than N66, they 
do contain high-excitation gas. 

4.2. PAH vs. warm dust 

We plot in Fig. [3] the PAH/14/im ratio as a function 
of the ISRF hardness in the three giant H ii regions of 
our sample. As explained in detail in the following, we 
distinguish 2 behaviors depending on the ISRF hardness. 
For a soft radiation field, the PAH/14 /im ratio is fairly 
independent on the IRSF hardness and shows a disper- 
sion of less than a factor 5. For a hard radiation field, 
PAII/14/xm shows a much wider dispersion with values 
that are systematically equal to or lower than the average 
ratio observed in a soft radiation field environment. 

Although not shown, we obtain a fairly similar trend 
when analyzing the sum of the PAH equivalent widths 
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TABLE 3 

Molecular hydrogen line fluxes. 





Ha S(l) 


Ha S(2) 


NGC3603#3 


1.62 (0.60) 


1.73 (0.60) 


NGC3603#4 


<1.57 


1.43 (0.60) 


NGC3603#5 


0.83 (0.60) 


1.07 (0.60) 


NGC3603#6 


<0.38 


3.49 (0.60) 


NGC3603#7 


<6.28 


3.53 (0.60) 


NGC3603#8 


<11.4 


<7.25 


30DOR#2 


0.37 (0.10) 


0.41 (0.10) 


30DOR#3 


0.12 (0.10) 


< 0.11 


30DOR#4 


<0.20 


0.18 (0.10) 


30DOR#5 


0.52 (0.10) 


0.46 (0.10) 


30DOR#6 


0.50 (0.10) 


0.44 (0.10) 


30DOR#7 


0.31 (0.10) 


0.22 (0.10) 


30DOR#8 


0.65 (0.10) 


0.60 (0.10) 


30DOR#10 


0.42 (0.10) 


0.53 (0.10) 


30DOR#11 


0.36 (0.10) 


0.41 (0.10) 


30DOR#12 


0.61 (0.10) 


0.59 (0.10) 


30DOR#13 


<0.10 


0.10 (0.10) 


30DOR#14 


0.68 (0.10) 


0.60 (0.10) 


30DOR#15 


0.80 (0.10) 


0.78 (0.10) 


30DOR#16 


<0.40 


<0.35 


30DOR#17 


0.50 (0.10) 


0.45 (0.10) 








N66#2 


0.07 (0.02) 


0.05 (0.02) 


N66#3 


0.05 (0.02) 


0.02 (0.02) 


N66#5 


0.09 (0.02) 


0.08 (0.02) 


N66#6 


0.15 (0.02) 


0.12 (0.02) 


N66#7 


0.15 (0.02) 


0.12 (0.02) 


N66#8 


0.14 (0.02) 


0.15 (0.02) 


N66#9 


0.29 (0.02) 


0.29 (0.02) 


N66#10 


0.29 (0.02) 


0.29 (0.02) 


N66#ll 


0.29 (0.02) 


0.29 (0.02) 


N66#12 


<0.07 


<0.10 


N66#13 


0.10 (0.02) 


0.13 (0.02) 


Note. — Line fluxes 


are expressed in > 


dO"^" Wcm"^. Tlio uncer- 



tainty is given in parentheses. 



§ 100 



+ 10 



















- NGC3503 — E 




- 30Dorodus ' 




- N66 





[Nelll] / [Nell] 

Fig. 2. — The quantity ([Ne ll]+C[Ne lll])/Huo (which is propor- 
tional to the neon abundance, see text) is plotted against the gas 
excitation as traced by the neon line ratio [Ne ni]/[Ne ll]. Squares 
indicate ionizing clusters, circles regions with silicate in emission, 
filled diamonds PDRs, open-diamonds ionized gas, and half- filled 
diamonds regions showing both PDR and ionized gas tracers (see 
Sect. 13.21 and Lebouteiller et al. 2008 for the spectra and the main 
spectral characteristics). 




1 10 
[Ne'll] / [Nell] 



Fig. 3. — The PAH emission measured against the lAfim dust 
continu um is plotted as a function of the ISRF hardness as defined 
in Sect. l3.2T3l The horizontal gray band illustrates the location of 
positions with significant PAH emission, while the vertical band 
indicates the threshold above which PAHs are believed to be de- 
stroyed by energetic photons. See Fig. [2] for the symbol explana- 
tion. 

rather than the PAH/14 /xm flux ratio. This is because 
the dust continuum lying under the PAH features is on 
first order correlated with the continuum at 14 /ini. By 
normahzing the strength of the PAH features with the 
14 /im continuum, we aim at minimizing systematic ef- 
fects due to the fact that the dust continuum between 
6 /im and 14 fim may show variations (different dust tem- 
perature, presence of silicate dust emission/absorption, 
stochastic heating of dust grains vs. thermal equihb- 
rium). 

4.2.1. Stellar clusters 

The central steflar clusters (position #6 in NGC 3603, 
#3 in 30 Dor, and #6 in N 66) lie somewhat off the trend 
identified in Sect. 14.21 They display low PAH/14 /im ra- 
tios — which is expected since PAHs are likely destroyed 
by photons from massive stars in the cluster — but the 
ISRF hardness is not the highest in the sample. This 
was also observed in NGC 3603 at small spatial-scales by 
Lebouteiller et al. (2007) who concluded that the pres- 
ence of a cavity around the central stellar cluster (volume 
in which [Ne ill] and [S iv] would strongly emit) cleared 
out by massive stars can explain this effect. In the case 
of N66, the most massive stars are actually not in the 
cluster NGC 346 and the ISRF hardness peaks to the east 
of it (Whelan et al. 2010). 

4.2.2. Soft radiation field environment 

The spectra toward all the positions characterized by 
soft radiation ([Ne iii]/[Ne ii]< 3) are PAH-dominated 
and correspond to individual PDRs (filled diamonds in 
Fig. |3|). The dispersion of the PAH/l4/xm ratio along 
the "PDR stripe" (around a factor 5, indicated by the 
horizontal band) is much smaller than the global varia- 
tions between all the positions (« 3dex), implying that 
the PAHs and the dust particles responsible for the con- 
tinuum emission scale with each other on first order in 
PDRs. For a soft radiation field, one expects the dust 
continuum in PDRs to be dominated by stochastically 
heated grains (very small grains; VSGs) (see e.g., Ra- 
pacioli et al. 2006; Berne et al. 2007). Given the size 
of the extraction window of the closest giant H ii region 



8 



Lebouteiller, Bernard-Salas, Whelan, Brandl, Galliano, Charmandaris, Madden, Kunth 



NGC3603 (Sect.|3II]), this suggests that PAHs and VSGs 
emit in a zone of less than 0.4 pc across (projected). 

Sources with significant 14 fim emission unrelated to 
the present star-formation episode are expected to lie 
well below the PDR stripe (horizontal band in Fig.|3]), 
while sources with significant PAH contamination from 
excitation by older stellar populations are expected to lie 
above. The upper limit on the potential contamination 
is given by the dispersion along the PDR stripe, i.e., a 
factor 5. Since the molecule and dust excitation in giant 
H II regions is dominated by massive stars, we consider 
the PDR stripe as a reference for interpreting integrated 
emission of galaxies where other sources of excitation can 
contribute (Sect. [5]). 

4.2.3. Hard radiation field environment 

Sources with low PAH/14/im ratios are systemati- 
cally characterized by hard radiation in the ionized gas 
([Ne iii]/[Ne ii]> 3; Fig.[3]). The PAH and the dust emis- 
sion thus become decoupled in extreme environments. 
Most of the positions with low PAH/14/im ratios were 
identified as pure ionized gas close to the central clusters 
by Lebouteiller et al. (2008) using both MIR imagery and 
spectroscopy (empty and half- filled diamonds in Fig. [3]). 
In the ionized gas, the dust continuum at 14 /um is ex- 
pected to be dominated by big grains in thermal equilib- 
rium with the radiation field. 

For a given giant H ii region, observed variations in 
the PAH/14/im ratio are not controlled by metallicity 
effects since metal abundances are found to be remark- 
ably uniform (Lebouteiller et al. 2008). The drop seen 
around [Ne lll]/[Ne ii]> 3 is due to a combination of two 
effects: the transition from VSGs in the PDRs to big 
grains in the ionized gas which results in increasing the 
lAfim continuum level (see also Dale et al. 2000), and 
the transition from PAH excitation in the PDRs to PAH 
destruction in the ionized gas. Although both effects are 
degenerate in Figure [3j a variation of the intrinsic PAH 
abundance is necessary to explain the observed variations 
traced by the PAH/14 ^m ratio (see Sect.|43]and Fig.[6]). 
Destruction of PAHs is most likely responsible for the 
low PAH abundance in the ionized gas, as it has already 
been observed in Galactic H ii regions (e.g., Geballe et al. 
1989; Verstraete et al. 1996; Cesarsky et al. 1996; Mad- 
den et al. 2006). Although it is also possible that the 
density profile controls the PAH abundance if the latter 
results from an equilibrium between thermal evaporation 
and reconstruction by encounters with (Giard et al. 
1994), our results suggest that the ISRF hardness plays 
an important role in the molecule destruction (see also 
the comparison between PAH and H2 emission in Sect.|6]). 

4.3. Dust emission properties in N 66 

The PAH/14/im values in N66 are larger than what 
could be expected from the radiation field properties, 
since all the sources in this object have [Ne iii]/[Ne 11] > 4 
(Fig. [3|). Figure|4] shows that continuum- normalized 
spectra of typical PDR-dominated sources in N66 (#10) 
and 30 Dor (#14, #15) have similar PAH integrated 
strength. It must be reminded here that the aperture 
physical size is almost identical for the observations in 
N66 and 30 Dor as these objects lie almost at the same 
distance fSect. l3TT|) . Although the PAH band ratios dif- 
fer in N66 and 30 Dor (see Galliano et al. 2008a), the 
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Fig. 4. — The positions #14 and #15 in 30 Dor (combined - 
green) show the closest continuum-normalized MIR spectrum to 
that of the position #10 in N 66 (blue). The continuum was re- 
moved using a spline function. The 11.3 fim PAH band across N 66 
is consistently stronger relative to the other bands as compared to 
the other H II regions. 

total PAH luminosity is identical. Hence the fact that 
PAH/14/xm values in N66 are offset with respect to 
30 Dor sources is due to relatively reduced warm dust 
emission in N66 rather than enhanced PAH emission. 
Could the weak dust emission in N 66 be a direct effect 
of metallicity on the dust grain abundance? The dust 
emission in the ionized gas is indeed expected to be re- 
duced in low-metallicity environments, since the dust-to- 
gas ratio scales directly with the ISM metallicity (Leroy 
et al. 2007; Lisenfeld & Ferrara 1998). If PAH emission 
would follow the same metallicity-dependence as dust 
emission, one would expect the PAH/14 /im ratio in N66 
and 30 Dor to be identical for a given ISRF hardness. 
Our results therefore seem to imply that PAH emission 
and dust emission do not scale together with metallicity 
and, more to the point, that PAH emission would de- 
pend less on metallicity than the warm dust emission. 
This is however difficult to reconcile with models of dust 
formation (Sect. lO]) . 

The relatively low metallicity of N66 could have an 
indirect influence via a change in the physical condi- 
tions leading to enhanced molecule photodestruction. 
Metal-poor sources are usually characterized by a rel- 
atively harder ISRF (Sect. [521) ^^'^ ^-l^o have a more in- 
tense ISRF because of the relatively low dust-to-gas ratio 
which reduces the ISM opacity. Nevertheless, NGC 3603 
and 30 Dor, which are more metal-rich than N66, also 
contain regions with a comparable excitation degree 
(Sect. 1322]). Therefore, if metallicity had only an indi- 
rect influence through variations of the ISRF properties, 
data points of the H 11 regions should align with each 
other in Fig.O 

Besides metallicity, a shift in the PAH/14 /im ratio can 
be observed because of different ISRF properties affect- 
ing the dust heating mechanism. For a given ISRF hard- 
ness ([Ne iii]/[Ne 11]), a lower dust temperature implies 
a higher PAH/14/im fSect. lO| . 

4.4. Theoretical influence of the ISRF properties on the 
PAH/14 fJ-m ratio 



Influence of the Environment on PAH Emission in Star-Forming Regions 



9 



Independently from metallicity and the possible con- 
tribution of different phases within a given beam, the 
simple variations of the illumination conditions of the 
PAH and dust grains can cause significant changes in 
the observed PAH/14 /im ratio. In order to quantify 
these variations, we have modeled realistic dust spec- 
tra with a homogeneous ISRF. We assumed the dust size 
distributions and abundances of the diffuse Galactic ISM 
modeled by Zubko et al. (2004), assuming bare grains 
with solar abundance constraints. This grain mixture is 
made of PAH (neutral and cationic; Draine & Li 2007), 
graphite (Laor & Draine 1993), and silicate (Weingart- 
ner & Draine 2001). The temperature fluctuations are 
computed using the method of Guhathakurta & Draine 
(1989). 

In order to parametrize the ISRF, Jx, we synthesize 
a single instantaneous burst with the stellar population 
synthesis model PEGASE (Fioc & Rocca-Volmerange 
1997). We measure the synthesized ISRF at Myr. Since 
PAHs emit principally from neutral regions (PDRs and 
diffuse ISM), we cut out the ionizing continuum of this 
ISRF. The family of ISRFs is parametrized by two quan- 
tities, its energy density: 

U= / Jxd\ (2) 

~'0.0912 /jin 

and its hardness: 

/ JxdX 

(H - Jr'''- • (3) 

/ Jx/{hc/X)d\ 

Jo. 0912 

The energy density is a simple scaling factor of the 
ISRF. The hardness parameter quantifies the weight of 
the UV photons compared to the redder ones. There 
are many ways to modify the hardness of an ISRF. We 
chose to apply a screen extinction exp(— t(A)), varying 
the optical depth. This choice is purely formal. High 
optical depths make the ISRF redder and decrease the 
subsequent {hv) parameter. 

The dust spectra generated are shown in Fig. [5] The 
top panel of Fig.[S] demonstrates the effect of the ISRF 
hardness on the shape of the dust SED. Hard ISRFs have 
prominent PAH features. When the hardness decreases 
the PAH strength decreases, as the UV-to-visible ratio 
decreases. This is due to the fact that PAHs have a 
much steeper cross-section in the UV, compared to other 
grain species (e.g.. Fig. 2 of GaUiano et al. 2008b). 

On the other hand, the bottom panel of Fig. [5] shows 
the variations of the shape of the dust spectrum with 
the energy density. In this case the hardness of the ISRF 
is fixed. PAHs are stochastically heated. They are con- 
trolled by single photon heating. It means that the shape 
of their emission spectrum will depend only on the av- 
erage energy of the incident photons {{hv)), and not on 
their rate {U). The very small grains, responsible for the 
flat continuum between 20 /im and 40 fim , for W < 10 
on the bottom panel of Fig.[5l are in the same regime. 
However, when the energy density increases, the equilib- 
rium temperature of the large grains increases, and the 
largest small grains reach thermal equilibrium. Conse- 
quently the emission peak shifts at shorter wavelengths. 




1 10 100 1000 

Wavelength [yum] 

Fig. 5. — Dust spectra used to study the effects of the illumina- 
tion conditions on the PAH/14 /xm ratio, from a purely theoretical 
point of view. Both panel represent families of dust spectra exposed 
to a single ISRF. Each spectrum is normalized by its bolometric 
luminosity. The incident ISRF is parametrized by its total energy 
density U and by its hardness. The top panel demonstrates the 
effect of the hardness, fixing U. The bottom panel demonstrates 
the effect of the starlight intensity, fixing the hardness. 

For high energy densities {hi > 10^), the grains respon- 
sible for the 14 /im emission are mainly at thermal equi- 
librium. Therefore, their emissivit\rn is enhanced, and 
the PAH/14/im ratio decreases. 

We refer to the monochromatic emissivity t^{u) (per unit 
frequency, per unit mass, averaged over solid angle) of dust 
grains at thermal equilibrium with the ISRF. In this context, 
= X By{Teq,v), where k,{v) is the grain opacity and 

the Bv(Tp^q,u) is the Planck function for an equilibrium tempera- 
ture Teg. 
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Fig. 6. — Theoretical variations of the PAH/14 /.im ratio, as a 
function of the illumination conditions. Those curves are a directly 
integrated from the SEDs of Fig. [5] U is the energy density of 
the ISRF, and (hu) its hardness. The top axis gives the effective 
temperature of a single O or B star that would have the hardness 
of the bottom axis. 

FigurelH shows the effect of varying (hv) and U on 
the PAH/14 /im ratio. For comparison with our mea- 
surements, the flux of the PAHs, /(PAH), includes all 
the bands between 5 /zm and 12 /xm. The 14 /im flux, 
/(14 /xm), is the dust continuum integrated between 
13.8 fim and 14.8 /im. We have indicated, on the top 
axis, the hardness corresponding to a single star (black 
body), with an effective temperature Tgtar- Figure|6] 
demonstrates that, in the range corresponding to the re- 
gions studied in this paper (Tatar > 10** K, or (hv) > 2; 
Fig.©, the spread in the PAH/14 /xm ratio due to vari- 
ations in the illumination conditions is at most a factor 
of « 2.6. It is therefore much lower than the observed 
spread 100), which consequently requires variations 
of the intrinsic abundance of PAHs. Our results suggest 
that PAH destruction in the ionized gas could explain 
the observed variations of the PAH/14/im ratio across 
each region (Sect. 14. 23)) . 

4.5. PAHs vs. ionized gas 

We next compare PAH emission to the H i recombi- 
nation line Hua which traces the warm photoionized gas 
(Sect. l33?2|) . While the 14 /im continuum can be domi- 
nated by a different dust component in PDRs and in the 
ionized gas, Hua simply provides a tracer of the ionized 
gas. As a result, the PAH/Hua ratio is thus propor- 
tional to the relative contribution of PDR and ionized 
gas within the extracted beam. 

Although the quantity PAH/Hua is a priori expected 
to depend on metallicity to first order because PAH 
abundance itself scales with carbon abundance, our re- 
sults show that the 14 /im dust continuum seems to 
be depend more on metallicity than the PAH emission 
fSect. l42|) . We therefore choose not to normalize the 
PAH/Hua ratio by the metallicity. Further discussion 




1 10 
[NePi] / [Ne"] 

Fig. 7. — The PAH emission (measured against Hua) is plotted as 
a function of the ISRF hardness. See Fig.[3]for symbol explanation. 

on the influence of metallicity on the dust continuum 
emission can be found in Sect. 15. 31 

The PAH/Hua ratio shows a smooth trend as a func- 
tion of the ISRF hardness (Fig.!?)). The main difference 
between the PAH/Hua ratio and the PAH/14 /im ratio 
(Fig. [3]) concerns the uniformity of the trend for the three 
giant H ii regions as a function of the ISRF hardness 
using PAH/Hua. In particular, the sources in N66 now 
cluster with those in NGC 3603 and 30 Dor, in agreement 
with Sect. l4.3l in that for the same physical conditions the 
14 /im dust continuum in this object is fainter relatively 
to the other giant H ii regions. This result implies that 
the three giant H ii regions follow the same transition 
from PDR dominated sources (filled diamonds in Fig. [71) 
to ionized gas dominated sources (empty diamonds) re- 
gardless of the environment metallicity. 

Since in star-forming regions PAHs and Hua trace re- 
spectively the PDRs and the ionized gas, we expect the 
PAH/Hua to reflect the mixture PDR/ionized gas. For a 
soft radiation field, the PAH/Hua ratio seems to reach a 
constant value. We attribute this effect to the size of the 
extraction window which makes it impossible to probe 
"pure" exposed PDRs with no ionized gas counterparts 
(such as in NGC 7023; Cesarsky et al. 1996), even in the 
closest object of our sample, NGC 3603 fSect. lO]) . The 
examination of individual spectra confirms that none of 
the regions resembles an exposed PDR (Lebouteiller et 
al. 2008). Although it is clear that ionized gas always 
mixes along the lines of sight from the observations of 
lines such as Hua, [Ne ii], or even [Ne ill], it is difficult 
to assert what fraction of it is physically associated with 
the PDRs probed. Nevertheless, we notice that PAH 
dominated regions in 30 Dor do not show systematically 
larger PAH/Hua ratios than in NGC 3603 despite the 
larger physical aperture, indicating that contamination 
by ionized gas in the beam is not significant (or that the 
contaminations by surrounding regions to the PAH and 
Hua emission are identical). 

4.6. Influence of shocks 

Shock conditions might modify the PAH abundance 
locally (Micelotta et al. 2010). O'Halloran et al. (2006) 
suggested that PAHs in low-metallicity dwarf galaxies 
could be destroyed by SN blast waves, based on the use 
of the [Fe ii]/[Ne ii] ratio. In shocked regions, the iron 
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Fig. 8. — The PAH/Hua ratio is plotted as a function of the 
abundance ratio Fe/Ne. See Fig. [2] for symbol explanation. 

locked onto dust grains is released into the ISM, increas- 
ing the iron abundance. As noticed by O'Halloran et al. 
(2006), large [Fe ii]/[Ne ii] ratios in low-metallicity envi- 
ronments can also be due to the fact that the amount of 
iron trapped in dust grains depends on the dust-to-gas 
ratio which itself depends on metallicity (e.g.; Lisenfeld 
& Ferrara 1998). 

In order to explore whether shocks can modify the PAH 
emission in giant H ii regions, we use the Fe/Ne abun- 
dance ratio instead of [Fe ii]/[Ne il] to avoid possible bi- 
ases from ionization effects, although the iron abundance 
determinations also bear significant errors driven by the 
ionization correction factors. Abundances are taken from 
Lebouteiller et al. (2008). In FigureH we plot PAH/Hua 
against the Fe/Ne ratio. It can be seen that there is 
no clear trend of decreasing PAH/ 14 /im with increasing 
Fe/Ne. 

The Fe/Ne ratio toward position #17 (corresponding 
to a young stellar object under the influence of nearby 
massive stars and a supernova) is much larger than the 
values toward the other sources in this region. Never- 
theless, the corresponding PAH/14/xm ratio toward this 
object is among the highest in our sample, indicating 
that even though shocks are present, they do not re- 
sult in significant PAH destruction. This result agrees 
with Indebetouw et al. (2009) who found that the ion- 
ization structure of 30 Dor is dominated by photoioniza- 
tion rather than shocks. The wide range of Fe/Ne in this 
object is likely driven by uncertainties on the abundance 
determination rather than by different shock conditions 
(Lebouteiller et al. 2008). 

We would like to stress that our results suffer from 
the fact that the PAH emission might not be asso- 
ciated with the shock itself (if present). For this 
reason, it seems also difficult to use MIR line ratios 
which are known to provide good constraints on shock 
conditions, such as [O iv]25.89/^m/[Ne ii]12.81/im, or 
[Fe ii]25.98/^m/[0 w]25.98 fim (Allen et al. 2008). 

5. STAR-FORMING GALAXIES 

The PAH equivalent width in star-forming galaxies has 
been shown to decrease in the presence of hard ISRF, us- 
ing ISO (Madden 2000; Madden et al. 2006; Thuan et al. 
1999) and Spitzer (Houck et al. 2004; Engelbracht et al. 
2005; Wu et al. 2006). In this section, we examine consis- 
tent measurements of a sample of star-forming galaxies 




Fig. 9. — The PAH/14 flux ratios is plotted as a function of 
the ISRF hardness as traced by [Ne lll]/[Ne III. The plot for the 
giant H II regions alone is presented in Figure[3l 

to compare to the results we obtained for the giant H ii 
regions. 

5.1. PAH vs. warm dust 

We plot in Figure IH] the relation between PAH/14/im 
and the ISRF hardness for starburst galaxies, H ii galax- 
ies, and BCDs together with those of the giant H ii re- 
gions. Galaxy samples are presented in Sect. 13. l1 Star- 
burst galaxies and H ii galaxies are characterized by 
relatively soft radiation, [Ne iii]/[Ne ii]< 3, which is 
expected given their high metallicity (Sect. [121) ■ In 
contrast, BCDs, as they are more metal-poor than 
starburst galaxies, have relatively harder radiation 
(1.5 <[Ne iii]/[Ne ii]< 18). 

Overall, the PAH/14 /im in starburst galaxies and H ii 
galaxies is equal or higher than values in PDRs of giant 
H II regions, being consistent with most of their MIR 
brightness originating from PDRs. Most of the disper- 
sion is likely driven by contributions from other galaxy 
components (ISM, embedded compact H ii regions, ex- 
posed PDRs, and refiection nebulae) to the PAH and 
dust emission (see also Peelers et al. 2004). 

The PAH/14 /im in BCDs spans all the range from 
PDR values to ionized gas values (Sect. 2]). Four BCDs 
have significantly low PAH/14 /im ratios (Mark 1450, 
NGC 1569, IIZw40, and Haroll), and they are all char- 
acterized by [Ne iii]/[Ne ii]> 3. In the case of Haroll, 
the remarkably low PAH/14/im ratio is likely related to 
the enhanced dust emission observed around 24 /im (En- 
gelbracht et al. 2008; Galametz et al. 2009). 

5.2. Mixture of PDR - ionized gas 

The PAH/Hua correlation with the ISRF hardness in 
star-forming galaxies provides a slightly different picture 
(Fig. [TOk). The outlier Haroll in Figure [9] now lines 
up with the other objects like it was observed for N66 
confirming peculiar warm dust emission in both objects 
fSect.l431). 

Starburst galaxies seem to have systematically larger 
PAH/Hua values than PDR sources in giant H ii re- 
gions. We confirm this result by using [Ne ii]-t-C[Ne ill] 
as a tracer of the ionized gas (Fig.llOb). The 
PAH/([Ne ii]-fC[Ne iii]) ratio should be less depe ndent 
on metallicity than PAH/Hua (see also Sect. 14. 5| ). and 
it allows us to increase the statistics by adding the H ii 
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Fig. 10.— The PAH/Hua (a) and PAH/([Ne ii]+C[Ne iii]) (b) 
flux ratios are plotted as a function of the ISRF hardness as traced 
by [Ne lll]/[Ne ll], 

galaxies for which no Hua measurement is available. 

These results suggest that a fraction of the PAH emis- 
sion in starburst galaxies does not originate from a col- 
lection of PDRs similar to the ones observed in giant H ii 
regions. The contribution from other galaxy components 
probably affects significantly the PAH emission because 
of the large physical aperture size. Peeters et al. (2004) 
found that the ratio between PAH luminosity and other 
massive star-formation tracers (Lyman continuum pho- 
tons, infrared luminosity) could vary by as much as a 
factor 10 between H II regions and integrated galaxies, 
which is on the same order of magnitude as the the shift 
in PAH/Hua we observe between starbursts and individ- 
ual PDRs in giant H ii regions. 

Most BCDs show PAH/Hua ratios comparable to the 
values found in the ionized gas of the H ii regions. The 
morphological resemblance between MIR spectra of ion- 
ized gas in H ii regions and spectra of BCDs has been 
pointed out by several studies (e.g., Galliano et al. 2003; 
2005; Madden et al. 2006; Peeters et al. 2004). Our re- 
sults are consistent with the conditions for PAH emission 
and destruction in galaxies being on first order similar to 
those seen in giant H ii regions. 

5.3. Warm dust emission properties 

Figure[Tl] shows the ratio between the 14 fj,m contin- 
uum and ionized gas tracers as a function of metallic- 
ity (Sect.[323)- Two ionized gas tracers are considered, 
[Ne ii]-|-C [Ne iii] and the H i recombination line Hua, 
with the latter having the advantage of not scaling di- 
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Fig. 11. — The 14 /xm dust grain emission is compared to 
[Ne ll]+C[Ne III] (a) and Huo (b) as a function of the metal- 
licity. For clarity, we only display the average ratio in each giant 
H II region, with the associated standard deviation. For the star- 
burst galaxies and BCDs, we display the metallicity (in solar units) 
derived from MIR lines with a symbol, and the difference to the 
optical determination (when it exists) is shown with a horizontal 
bar. The black solid line shows the polynomial fit to all the data 
points. 



rectly with the abundance of heavy elements (Sect. [3T2^ . 
We assume in the following that most of the 14 dust 
emission originates in the ionized gas (Sect. lSmi) . 

Contrary to PAH/14 /im or PAH/Hua, the ratio be- 
tween the 14 /im continuum and ionized gas tracers is 
fairly insensitive to the ISRF hardness. Considering the 
combined results of the giant H il regions and the star- 
forming galaxies, the 14/xm/([Ne ii]-|-C [Ne ill]) ratio 
does not show any significant correlation with metallic- 
ity (Fig.inji). 

This contrasts with the 14/im/Hua ratio, especially 
for low-metallicity regions (Fig.lTTb). We conclude that 
the rather constant observed 14/im/([Ne ii]-|-C [Ne III]) 
ratio is due to the fact that the dust continuum scales 
with metallicity. Unfortunately, the large dispersion 
hampers a quantitative determination of the metallicity- 
dependence. 

Our results showed that the PAH/14 /im ratio in N66 
is relatively higher when compared to positions with sim- 
ilar ISRF hardness in other H Ii regions (Sect. 2]). This 
discrepancy is consistent with the warm dust continuum 
emission being reduced in this object because of the rela- 
tively low-metallicity environment. The dependence can 
be explained by a lower dust mass at low-metallicity (see 
Engelbracht et al. 2008). 
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6. PAHS VS. MOLECULAR HYDROGEN 

Emission from PAH molecules and molecular hydro- 
gen both originates close to or in PDRs (Sect l2.2( ). Since 
they are carbon-based molecules, the total abundance of 
PAHs in a galaxy is expected to depend at least partly on 
the metallicity. On the other hand, H2 formation on dust 
grains also implies an indirect metallicity-dependence. 
PAH and H2 survival in harsh radiation field bears im- 
portant differences since photodissociation of H2 requires 
energies >11.2eV, while much higher energies are needed 
for PAHs which are able to thermalize rapidly due to ef- 
ficient electron delocalization (Sect. [2^ . 

H2 transitions in the MIR give access to the bulk of 
the warm molecular phase (few 100 K). The most promi- 
nent MIR lines are due to the spontaneous deexcita- 
tion toward lower rotational levels within the ground 
vibrational-electronic state. Excitation of H2 can oc- 
cur in several ways, among which the most frequent 
are pumping of non-ionizing photons followed by fluores- 
cence, coUisional excitation of the lower rotational levels 
for large densities (few 10'^ cm~'^), thermal excitation via 
collisions with photoelectrons, and heating by shocks (see 
e.g., Roussel et al. 2007). 

In UV-dominated regions, PAHs and H2 are expected 
to trace each other at sufficiently large spatial-scales. 
Roussel et al. (2007) studied H 11 nuclei within the 
SINGS sample (Kennicutt et al. 2003) and found a 
tight correlation between PAH emission (bands within 
the Spitzer/IRAC [8/ini]band) and H2 emission (S(0)-|- 
S(l)-hS(2) lines). The PAH/H2 ratio has been found by 
these authors to be fairly insensitive to the radiation field 
intensity, which reinforces the hypothesis of a similar ori- 
gin in PDRs. Although with a smaller range of radiation 
field intensity, Bernard-Salas et al. (2009) obtained the 
same result in their sample of starburst galaxies. 

In order to trace the PAH emission in the objects of 
our sample, we use the most prominent bands (6.2 /im, 
7.7/im, 8.6 //m, and 11.3 /im), while we use the S(l) and 
S(2) lines to trace H2 (Sect. [3XT|) . Figure [12] shows 
the resulting PAH/H2 ratio as a function of the ISRF 
hardness. On first order, the PAH and H2 emission scale 
together if one considers that most of the positions with 
low PAH/H2 ratios are in fact lower limits. Only a small 
fraction of the dispersion (which is about a few tens) 
should be due to the lack of detection of the S(0) line, 
contributing about 25% to the total of S(0-2) in H 11 
nuclei fSect. l3.2TTl Roussel et al. 2007). A more detailed 
analysis reveals that the values in NGC 3603 and 30 Dor 
lie at the lower end of the range observed for the starburst 
galaxies, while the BCD Mark 1450 and most positions 
in N66 are characterized by significantly lower PAH/H2 
ratios. 

The PAH/H2 ratio shows a factor 10 variation among 
the starburst galaxies. This large dispersion can be at- 
tributed to contribution of older stellar populatons to the 
PAH emission, agreeing with results in Sect. 15. 21 

We notice that low PAH/H2 values are only found in 
objects with hard ISRF, reminding the results obtained 
in Sect.m Hence, it seems that PAH molecules are de- 
stroyed in these objects while H2 is still able to emit. 
Such a result is consistent with H2 being self-shielded so 
that it is able to emit regardless of the UV field intensity 
or hardness, contrary to PAH emission which is highly 
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Fig. 12. — The PAH emission is compared to H2 as a function of 
the ISRF hardness. See Fig. [2] for the symbol explanation. 

dependent on the physical conditions in the ionized gas. 
Note that this result supposes that the spectral extrac- 
tion window (Sect. 13. ip is large enough to include the 
contribution of both PAH-dominated and H2 dominated 
regions. 

The weak correlation of PAH/H2 with the ISRF hard- 
ness is less obvious than what we observed by using 
PAH/14 /im or PAH/Hua since for a hardness above the 
threshold for PAH destruction (Sect.Hl), the PAH/H2 ra- 
tio appears to be quite unrelated to the ISRF. Therefore, 
although these results strongly suggest that PAH and H2 
molecules depend differently on the incident UV field, an- 
other parameter might explain the range of PAH/H2 for 
similarly hard ISRF values. 

7. SUMMARY AND CONCLUSIONS 

In this paper, we first investigated the PAH emission 
properties of several sources within three giant H 11 re- 
gions, NGC 3603 in the Milky Way, 30 Dor in the LMC, 
and N66 in the SMC. The sample spans an interesting 
range of metallicities (down to 0.2 Z0) and physical con- 
ditions (ISRF hardness). The sample of giant H 11 regions 
contain several ISM components contributing to the MIR 
spectra of star-forming galaxies, stellar clusters, individ- 
ual PDRs, warm photoionized g well as embedded 
MIR bright regions. The main results are: 

• The PAH/14/im ratio shows a strong dependence 
on the ISRF hardness as traced by [Ne iii]/[Ne 11]. 
The locations of the data points in the PAH/ 14 /im 
vs. [Ne iii]/[Ne 11] diagram agree well with the PDR 
vs. ionized gas nature of the sources. For soft ISRF, 
PAH/14 ^m is uniform and dominated by PDRs. A 
threshold occurs around [Ne iii]/[Ne 11] ~ 3 above 
which the dispersion of PAH/14/i,m values greatly 
increases, with a decreasing ratio as the ISRF hard- 
ens. 

• N 66 is characterized by large PAH/14 /im values 
compared to positions in 30 Dor with similar phys- 
ical conditions. The difference is due to suppressed 
warm dust emission in N66, which we propose is 
an effect of the low dust-to-gas ratio in the ionized 
gas. 

• The variations of the PAH/14/ini ratio are not 
dominated by varying ISRF conditions leading to 
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different heating mechanism of PAHs and of the 
dust. A change in the PAH abundance is reqiiired, 
being consistent with the paradigm of photode- 
struction of PAH molecules in the ionized gas. 

• PAH emission is compared to the H i recombina- 
tion line Hua which traces the ionized gas. The 
variations of PAH/Hua are directly related to the 
relative contribution of PDR vs. ionized gas. We 
find that all sources in the giant H ii regions follow 
the same trend as a function of the ISRF hardness. 
We use these results as a reference for discussing 
observations of star-forming galaxies. 

• By investigating the dependence of PAH/Hua on 

the iron abundance, we do not find any evidence 
for enhanced PAH destruction by shocks. 

The results on the giant H ii regions are compared to 
those of star-forming galaxies (BCDs, starburst galaxies, 
and H II galaxies), with the following findings: 

• The PAH/14/xm ratio in star- forming galaxies fol- 
lows on first order the trend of the giant H II 
regions. Starburst galaxies globally lie with the 
PDRs in giant H ii regions while BCDs lie with 
sources in giant H II regions dominated by ionized 
gas. 

• PAH/Hua in starburst galaxies is larger than PDR 
values in giant H II regions which is likely due to 

the PAH excitation by old stellar populations in 
the diffuse ISM. The PAH emission component in 
the galaxies of the sample is thus not dominated by 



a collection of PDRs similar to the ones observed 
in giant H II regions. 

• We find a metallicity dependence of the 14 /xm/Hua 
ratio, which reflects the warm dust abundance in 
the photoionized gas. 

Low-metallicity star-forming galaxies are particularly 
challenging environments for detecting PAH emission be- 
cause of the low carbon abundance which limits the for- 
mation of hydrocarbon carriers, and of the hard and 
intense radiation field which is able to destroy PAH 
molecules. Furthermore, dwarf galaxies do not contain 
significant amounts of GMCs, implying little integrated 
volume of PDRs. Our findings suggest that metal-poor 
star- forming galaxies such as BCDs have PAH/Hua ra- 
tios that agree on first order with what is expected from 
typical physical conditions similar to the ionized gas in 
H II regions. Even though PAHs could be formed in large 
amounts over the galaxy history, molecules are mostly 
photodissociated in the ionized gas. In more metal-rich 
galaxies such as starburst galaxies, there is a sign of sig- 
nificant diffuse PAH emission which could severely bias 
star-formation rate indicators based on PAH emission. 
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